41)713172 


AFCRL-7WJ460 
19  AUGUST  1970 

AIR  FORCE  SURVEYS  IM  GEOPHYSICS.  NO.  224 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 


L.  G.  HAN5COM  FIELD.  BEDFORD,  MASSACHUSETTS 


Earth  Environmental  Noise  Fields 


FA.  CROWLEY 
H.A.  OSSJNG 


\ 


\ 


United  States  Air  Force 


Qualified  raquattore  may  obtain  additional  copies  from  tfw 
Detente  Documentation  Cantar.  Ail  othere  rbouid  apply  to  the 
Clearinghouse  for  Fadaral  Scientific  and  Technical  Information. 


AFCRL-70-0460 
19  AUGUST  1970 

AIR  FORCE  SURVEYS  IN  GEOPHYSICS,  NO.  224 


V  ’’  -if* 

TERRESTRIAL  SCIENCES  LABORATORY  PROJECT  7639 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 

L.  G.  HANSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


Earth  Environmental  Noise  Fields 


F.A.  CROWLEY 
H.A.  OSSING 


This  document  has  been  approved  for  public  release 
and  sole;  its  distribution  is  unlimited* 


* 


United  States  Air  Force 


Abstract 


A  new  method  to  approach  earth  noise  field  problems  in  inertial  component 
testing  is  presented.  This  method  is  considered  superior  to  accepted  "quiet  site" 
and  "motion  suppression  structure"  techniques.  The  approach  is  to  describe 
earth  noise  by  dt  lermining:  (1)  the  class  of  the  random  motion;  (2)  the  level  of 
motion;  (3)  the  mode  of  motion;  and  (4)  the  component  relationships.  The 
efficacy  and  rationale  of  this  technique  is  demonstrated  in  an  experiment  that  re¬ 
sults  in  the  removal  of  earth  noise  error  terms  from  gyro  data. 
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Earth  Environmental  Noise  Fields 


It  is  now  widely  recognized  t  further  advancements  in  the  accuracy  of 
inertial  guidance  systems  are  conditioned  on  properly  resolving  between  errors 
due  to  accelerometer  and  gyro  deficiencies  and  those  caused  by  false  assumptions 
about  the  nature  of  the  test  environment.  No  longer  can  component  performance 
tests  treat  the  earth  as  a  rigid  body  spinning  on  a  fixed  axis.  Indeed,  as  we  look 
ahead,  we  find  that  the  earth  is  a  noisy  medium  for  inertial  component  testing. 

Historically,  the  problem  of  earth  noise  in  inertial  component  tests  has  been 
dominated  by  two  opposing  test  philosophies.  At  the  one  extreme  there  are  those 
who  would  wish  to  avoid  earth  noise  problems  by  seeking  out  stable,  low-motion 
sites  for  locating  test  facilities.  At  the  other  extreme,  there  are  those  who  would 
rely  solely  on  motion  attenuation/compensation  structures  to  obtain  a  permissible 
test  environment.  To  these  two  approaches  we  now  add  a  third,  namely  the  sup¬ 
pression  of  earth  noise  terms  through  optimum  processing  of  auxiliary  earth 
motion  measurements.  In  effect,  we  seek  to  treat  earth  noise  as  an  extension  of 
shake  table  testing.  We  propose  to  do  this  by  measuring  earth  noise  and  then 
allowing  for  it  in  component  tests. 

AFCRL  was  introduced  to  the  problem  of  earth  noise  in  inertial  components 
by  the  NASA  Electronics  Research  Center  (ERC),  formerly  located  in  Cambridge. 
It  was  NASA's  attitude  that  the  advantage  of  being  an  abutter  to  MIT  outweighed 
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the  difficulties  of  constructing  structures  for  attenuating/compensating  in  the  high 
noise  levels  at  its  Cambridge  site.  In  con*r,ast,  the  Air  Force  sought  to  forestall, 
at  least  temporarily,  earth  noise  problems  by  going  to  a  seismically  quiet  site  in 
New  Mexico, 

In  their  extremes,  both  test  philosophies  are  found  wanting.  Those  that 
sought  out  low  noise  environments  discovered  that  test  facilities  bring  with  them 
motion  sources  that  degrade  their  once  "quiet"  site.  In  turn,  advocates  of 
motion- suppressing  structures  have  met  with  only  partial  success  in  suppressing 
earth  noise  over  the  pass  band  of  interest. 

Earth  noise  in  component  testing  can  be  defined  as  those  unknown  or  un¬ 
accounted  motions  of  the  test  environment  that  affect  inertial  hardware  perform¬ 
ance.  The  total  motion  at  a  test  facility  can  be  conceived  as  the  composite 
motion  due  to  propagating  waves,  deformations  due  to  loading,  local  thermal 
stresses,  alterations  in  ground  water,  etc. ,  as  well  a°  to  large  scale  phenomena 
like  wobble  of  the  earth's  axis  and  tectonic  processes.  For  our  part,  we  have 
concentrated  our  efforts  on  only  a  part  of  the  problem— namely,  motions  between 
0.1  and  10.0  Hz.  This  choice  of  frequencies  constitutes  a  troublesome  cross¬ 
over  zone  between  the  effective  bands  of  passive  and  active  motion  reduction 
structures  conceived  by  ERC-NASA,  our  original  sponsor. 

The  motions  at  the  bulk  of  the  sites  that  we  have  studied  support  the  thesis 
that  earth  noise  in  the  band  0. 1  to  10.0  Hz  takes  the  form  of  propagating  waves 
that  are  excited  by  a  large  number  of  independent  random  motion  sources.  The 
features  of  these  motions  relative  to  testing  are  determined  by  the  behavior  and 
distribution  of  the  sources  as  well  as  by  the  ground  transmission  characteristics 
of  an  area.  Measured  distributions  of  the  motion  (Figure  1)  have  the  form  of  a 
normal  distribution.  Normally  distributed  variates  will  result  when  the  measure¬ 
ments  are  the  sum  of  a  large  number  of  independent  random  sources.  The  tem¬ 
poral  attributes  of  normal  motion  processes  are  well  described  in  terms  of 
correlation  and/or  spectral  estimates. 

Figure  2  is  a  set  of  spectra  formed  on  data  taken  at  the  ERC-NASA  site,  over 
a  period  of  one  week.  The  levels  of  measured  spectra  are  found  to  be  strongly 
locked  to  the  work  week  cycle  for  this  urban  area,  in  geological  cross  sect! 
the  site  consists  of  an  alluvial,  slow-speed  layer  overlying  a  hard,  high-spied 
rock.  For  such  layering,  we  can  expect  strong  ground  resonances  to  occur.  The 
observed  narrow  band  spectra  appear  to  be  an  attribute  of  man-made  sources  and 
resonant  structure  of  the  site.  When  man-made  motion  sources  drive  local 
ground  resonances,  large  motions  result.  The  motion  levels  found  at  this  site  are 
several  orders  of  magnitude  larger  than  those  measured  at  local,  surburban, 
hard-rock  areas.  Our  motion  statistics  can  be  ♦-eated  as  periodically  stationary 
for  this  locale;  that  is,  correlations  are  insensitive  to  discrete  shifts  in  their 
time  of  estimation. 
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Figure  1. 
Density 


Acceleration  Probability 


The  assessment  of  earth 
motion  in  inertial  component  tests 
calls  for  establishing  two  other 
properties  of  the  motion  environ¬ 
ment  in  addition  to  class  and  level, 
namely  the  mode  and  the  organiza¬ 
tion  of  the  motion.  By  mode,  it  is 
meant  that  we  must  differentiate 
between  rotational  and  translational 
motions.  In  turn,  estimates  of 
organization  call  for  the  strength 
of  amplitude /phase  relations  be¬ 
tween  components.  Let  us  first 
consider  the  problem  of  isolating 
rotations  and  translations.  The 
problem  is  relevant  to  systems 
using  single-degree-of-freedom 
(SDF)  gyros.  SDF  gyro-based 
Azimuth  Laying  Sets  (ALS)  are 
now  being  deployed  as  the  prime 
azimuth  control  for  the  Minuteman 
fleet. 


The  problem  of  isolating  rotations  and  translations  lies  in  the  fact  that  our 
measuring  instruments,  pendulous  seismometers,  respond  equally  well  to  both 
rotations  and  translations.  To  overcome  this  ambiguity  in  measurement,  we  re¬ 
quire  the  combination  of  four  seismometers  arranged  as  in  Figure  3.  We  must 
combine  two  difference  motions  as  given  in  Eq.  (3)  of  the  figure. 


Figure  2.  Spectrogram,  ERC-NASA  Site  -  Kendall  Square, 
Cambridge,  Maas.  February  8-15,  1967 
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Figure  3.  Seismometer  Configuration  for 
Rotation,  Strain  Measurements  nt  the 
Honeywell  Facility 
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Figure  4.  Total  anil  Rotational  Motion 
Estimates  for  the  Honeywell  Test 


In  principle,  this  technique 
presents  no  real  difficulty;  in 
p"actice,  it  requires  highly  matched 
transducers,  because  linear  mo¬ 
tions  typically  are  much  larger 
than  rotations.  Using  this  approach, 
we  estimated  rotation  levels  on 
stands  now  being  used  for  ALS 
acceptance  tests.  Transducer 
matching  was  accomplished  off¬ 
line  using  Wiener  filters.  The 
results  of  this  experiment  are  given 
in  Figure  4.  Here,  is  shown  the 
mean-square  value  of  the  linear 
and  rotational  motions  as  a  func¬ 
tion  of  bandwidth.  Our  approach 
appears  to  be  adequate  for  es¬ 
timating  rotations  when  rotations 
approach  bothersome  levels. 

Having  discussed  the  mode  of 
the  motion,  we  now  turn  to  esti¬ 
mates  of  component  organization. 
Error  terms  introduced  by  non- 
linearities  in  both  gyros  and  pen¬ 
dulous  accelerometers  are  affected 
by  component  motion  relations. 

For  the  pendulum,  the  error  is 
known  as  zero  shift;  for  the  gyro, 
this  class  of  error  is  called  coning 
or  kinematic  rectification.  Es¬ 


timates  of  coning  e  t  or  in  an  SDK  gy  ro  cal!  for  computing  the  cross -correlation 
between  the  components  of  stand  rotations  about  the  gyro's  output  and  input  axes. 

A  considerable  complication  occurs  in  making  coning  error  estimates  if  *he  stand 
motion  is  not  a  normal  process,  for  then  i.igher-order  correlations  and/or  spectra 
must  be  determined. 

The  relations  between  component  motions  are  conveniently  revealed  through 
coherency  estimates,  a  measure  of  the  strength  of  component  relations  under 
spectral  decomposition.  Our  studies  show  the  organization  of  the  motion  to  be 
highly  site-dependent  To  illustrate  (Figure  5),  we  find  component  relations  are 
weak  on  Honeywell's  Aids  test  stands.  Here,  the  motion  environment  exhibits  the 
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Figure  6.  Northrop  Motion  Spectra 
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Figure  7.  Inside-Outside  Motion  (Time 
Domain) 


attributes  of  an  isotropic  motion 
field,  that  is,  the  computed  rela¬ 
tions  are  the  same  as  those  that 
would  result  from  a  set  of  inde¬ 
pendent  random  motion  sources 
uniformly  dispersed  around  the 
measurement  point.  On  the  other 
hand,  similarly  treated  measure¬ 
ments,  taken  at  Northrop's  test 
facility,  Norwood,  Mass.,  reveal 
that  the  motions  at  this  test  faci¬ 
lity  are  highly  organized  (Figure  6). 
Indeed,  when  we  use  an  optimum 
least- squares  prediction  to  remove 
the  motions  caused  by  exterior 
sources,  we  find  that  motions 
excited  by  intra-plant  sources  are 
well  related  and  concentrated 
around  8  Hz  (Figure  7). 

To  summarize,  in  component 
testing  a  description  of  earth  noist 
over  the  band  0.1  to  10.  0  Hz  calls 
for  determining  the  following. 

(1)  the  class  of  the  random  motion; 

(2)  its  level;  (3)  the  mode  of  the 
motion;  (4)  component  relationships. 

We  would  now  like  to  turn  to 
an  experiment  that  removes  earth 
noise  error  terms  from  gyro  data. 

In  effect,  we  treat  the  gyro  as  a 
multi-input  single-output  black  box 
to  earth  noise  incuts.  The  rationale 
of  the  experiment  is  outlined  in 
Figure  8. 

Figure  &  gives  a  sample  rela¬ 
tion  between  the  gyro  output  and 
one  of  the  component  motions  The 
top  portion  of  the  figure  is  the 
measured  gyro  spectra.  The  dark¬ 
ened  area  at  the  base  of  this  spectra 
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^ — i  coherency  is  t^at  portion  of  the  gyro  data  that 

i - 1 - 1  can  be  removed  by  seismic  measure- 

FREQUENCY  (Hi! 

ments,  using  a  non-realizable  Wiener 
Figure  9.  Gyro  Output  Spectra  Based  operator.  The  second  element  of  the 

on  Auxiliary  Motion  Estimates  figure  is  the  coherency  between  the 

data  pair.  In  order  to  evaluate  the 
coherency  plot,  we  must  consider 

L  levels  of  significance.  Figure  10 

total  ^c^dicteci  gives  the  expected  upper  bound  co¬ 

herency  estimate  for  the  conditions 
of  our  experiment  vs  bandwidth  at  a 
1  percent  level  when  the  computed 

o  ~  5  coherencies  are  actually  founded  on 

\ - 1 _ l 

frequency  (Hz)  an  uncorrelated  data  pair.  The  mea¬ 

surements  shown  in  the  previous 

i  100  figure  (Figure  9)  are  barely  signifi- 

&  I  1  di  1 1  cant.  In  turn.  Figure  11  is  the  co- 

}  herency  computed  between  accelera- 

o  5  tions  squared  terms  and  the  gyro  out- 

a  < _ J 

frequency (hz )  put.  Here<  the  computed  coherencies 

are  significant  at  well  in  excess  of  a 

Prediction  1  Percent  level-  In  our  experiment, 

we  find  we  can  remove  something 
like  50  percent  of  the  low  frequency 

content  of  the  gyro  output.  Compliance  in  the  gyro  is  known  to  lead  to  errors 
sensitive  to  acceleration  squared  terms. 

Heartened  by  this  success,  we  have  initiated  a  comprehensive  long  term  ex¬ 
periment  to  develop  an  estimate  of  the  ultimate  sensitivity  of  the  SDF  gyro  used  in 
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Figure  9.  Gyro  Output  Spectra  Based 
on  Auxiliary  Motion  Estimates 
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the  ALS  system  when  all  environmental  noise  terms  are  removed.  The  experiment 
is  a  joint  effort  by  AFCRL  and  SAMSO.  It  will  embrace  all  earth  noise  terms  to 
the  extent  that  we  can  recognize  and  account  for  them,  rather  than  just  the  0.1  to 
10.0  Hz  motions,  just  discussed.  The  experiment,  it  is  hoped,  will  have  a  marked 
impact  on  the  performance  loss  that  too  often  occurs  when  we  take  a  system  out  of 
a  contractor's  laboratory  environment  into  the  field.  It  should  also  point  the  way 
for  better  resolving  between  those  errors  due  to  earth  noise  and  those  due  to  hard¬ 
ware  deficiencies  in  future  inertial  component  tests. 
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